Current trends in RNA virus detection through metatranscriptome sequencing data.

With advances in sequencing technology, metatranscriptome sequencing from a variety of environmental and biological sources has revealed the existence of various previously unknown RNA viruses. This review presents recent major RNA virome studies sampled from invertebrate and vertebrate species as well as aquatic environments. In particular, we focus on severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and related RNA virus identification through metatranscriptome sequencing analyses. Recently developed bioinformatics software and databases for RNA virus identification are introduced. A relationship between newly identified RNA viruses and endogenous viral elements in host genomes is also discussed.

Genomic diversity of SARS-CoV-2 can be accelerated by mutations in the nsp14 gene.

Coronaviruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), encode a proofreading exonuclease, nonstructural protein 14 (nsp14), that helps ensure replication competence at a low evolutionary rate compared with other RNA viruses. In the current pandemic, SARS-CoV-2 has accumulated diverse genomic mutations including in nsp14. Here, to clarify whether amino acid substitutions in nsp14 affect the genomic diversity and evolution of SARS-CoV-2, we searched for amino acid substitutions in nature that may interfere with nsp14 function. We found that viruses carrying a proline-to-leucine change at position 203 (P203L) have a high evolutionary rate and that a recombinant SARS-CoV-2 virus with the P203L mutation acquired more diverse genomic mutations than wild-type virus during its replication in hamsters. Our findings suggest that substitutions, such as P203L, in nsp14 may accelerate the genomic diversity of SARS-CoV-2, contributing to virus evolution during the pandemic.

The SARS-CoV-2 spike S375F mutation characterizes the Omicron BA.1 variant.

Recent studies have revealed the unique virological characteristics of Omicron, particularly those of its spike protein, such as less cleavage efficacy in cells, reduced ACE2 binding affinity, and poor fusogenicity. However, it remains unclear which mutation(s) determine these three virological characteristics of Omicron spike. Here, we show that these characteristics of the Omicron spike protein are determined by its receptor-binding domain. Of interest, molecular phylogenetic analysis revealed that acquisition of the spike S375F mutation was closely associated with the explosive spread of Omicron in the human population. We further elucidated that the F375 residue forms an interprotomer pi-pi interaction with the H505 residue of another protomer in the spike trimer, conferring the attenuated cleavage efficiency and fusogenicity of Omicron spike. Our data shed light on the evolutionary events underlying the emergence of Omicron at the molecular level.

Efficient compression of SARS-CoV-2 genome data using Nucleotide Archival Format.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome data are essential for epidemiology, vaccine development, and tracking emerging variants. Millions of SARS-CoV-2 genomes have been sequenced during the pandemic. However, downloading SARS-CoV-2 genomes from databases is slow and unreliable, largely due to suboptimal choice of compression method. We evaluated the available compressors and found that Nucleotide Archival Format (NAF) would provide a drastic improvement compared with current methods. For Global Initiative on Sharing Avian Flu Data's (GISAID) pre-compressed datasets, NAF would increase efficiency 52.2 times for gzip-compressed data and 3.7 times for xz-compressed data. For DNA DataBank of Japan (DDBJ), NAF would improve throughput 40 times for gzip-compressed data. For GenBank and European Nucleotide Archive (ENA), NAF would accelerate data distribution by a factor of 29.3 times compared with uncompressed FASTA. This article provides a tutorial for installing and using NAF. Offering a NAF download option in sequence databases would provide a significant saving of time, bandwidth, and disk space and accelerate biological and medical research worldwide.

The SARS-CoV-2 Lambda variant exhibits enhanced infectivity and immune resistance.

SARS-CoV-2 Lambda, a variant of interest, has spread in some South American countries; however, its virological features and evolutionary traits remain unclear. In this study, we use pseudoviruses and reveal that the spike protein of the Lambda variant is more infectious than that of other variants due to the T76I and L452Q mutations. The RSYLTPGD246-253N mutation, a unique 7-amino acid deletion in the N-terminal domain of the Lambda spike protein, is responsible for evasion from neutralizing antibodies and further augments antibody-mediated enhancement of infection. Although this mutation generates a nascent N-linked glycosylation site, the additional N-linked glycan is dispensable for the virological property conferred by this mutation. Since the Lambda variant has dominantly spread according to the increasing frequency of the isolates harboring the RSYLTPGD246-253N mutation, our data suggest that the RSYLTPGD246-253N mutation is closely associated with the substantial spread of the Lambda variant in South America.

Enhanced fusogenicity and pathogenicity of SARS-CoV-2 Delta P681R mutation.

During the current coronavirus disease 2019 (COVID-19) pandemic, a variety of mutations have accumulated in the viral genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and, at the time of writing, four variants of concern are considered to be potentially hazardous to human society1. The recently emerged B.1.617.2/Delta variant of concern is closely associated with the COVID-19 surge that occurred in India in the spring of 2021 (ref. 2). However, the virological properties of B.1.617.2/Delta remain unclear. Here we show that the B.1.617.2/Delta variant is highly fusogenic and notably more pathogenic than prototypic SARS-CoV-2 in infected hamsters. The P681R mutation in the spike protein, which is highly conserved in this lineage, facilitates cleavage of the spike protein and enhances viral fusogenicity. Moreover, we demonstrate that the P681R-bearing virus exhibits higher pathogenicity compared with its parental virus. Our data suggest that the P681R mutation is a hallmark of the virological phenotype of the B.1.617.2/Delta variant and is associated with enhanced pathogenicity.

SARS-CoV-2 spike L452R variant evades cellular immunity and increases infectivity.

Many SARS-CoV-2 variants with naturally acquired mutations have emerged. These mutations can affect viral properties such as infectivity and immune resistance. Although the sensitivity of naturally occurring SARS-CoV-2 variants to humoral immunity has been investigated, sensitivity to human leukocyte antigen (HLA)-restricted cellular immunity remains largely unexplored. Here, we demonstrate that two recently emerging mutations in the receptor-binding domain of the SARS-CoV-2 spike protein, L452R (in B.1.427/429 and B.1.617) and Y453F (in B.1.1.298), confer escape from HLA-A24-restricted cellular immunity. These mutations reinforce affinity toward the host entry receptor ACE2. Notably, the L452R mutation increases spike stability, viral infectivity, viral fusogenicity, and thereby promotes viral replication. These data suggest that HLA-restricted cellular immunity potentially affects the evolution of viral phenotypes and that a further threat of the SARS-CoV-2 pandemic is escape from cellular immunity.

Sarbecovirus ORF6 proteins hamper induction of interferon signaling.

The presence of an ORF6 gene distinguishes sarbecoviruses such as severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 from other betacoronaviruses. Here we show that ORF6 inhibits induction of innate immune signaling, including upregulation of type I interferon (IFN) upon viral infection as well as type I and III IFN signaling. Intriguingly, ORF6 proteins from SARS-CoV-2 lineages are more efficient antagonists of innate immunity than their orthologs from SARS-CoV lineages. Mutational analyses identified residues E46 and Q56 as important determinants of the antagonistic activity of SARS-CoV-2 ORF6. Moreover, we show that the anti-innate immune activity of ORF6 depends on its C-terminal region and that ORF6 inhibits nuclear translocation of IRF3. Finally, we identify naturally occurring frameshift/nonsense mutations that result in an inactivating truncation of ORF6 in approximately 0.2% of SARS-CoV-2 isolates. Our findings suggest that ORF6 contributes to the poor IFN activation observed in individuals with coronavirus disease 2019 (COVID-19).

SARS-CoV-2 ORF3b Is a Potent Interferon Antagonist Whose Activity Is Increased by a Naturally Occurring Elongation Variant.

One of the features distinguishing SARS-CoV-2 from its more pathogenic counterpart SARS-CoV is the presence of premature stop codons in its ORF3b gene. Here, we show that SARS-CoV-2 ORF3b is a potent interferon antagonist, suppressing the induction of type I interferon more efficiently than its SARS-CoV ortholog. Phylogenetic analyses and functional assays reveal that SARS-CoV-2-related viruses from bats and pangolins also encode truncated ORF3b gene products with strong anti-interferon activity. Furthermore, analyses of approximately 17,000 SARS-CoV-2 sequences identify a natural variant in which a longer ORF3b reading frame was reconstituted. This variant was isolated from two patients with severe disease and further increased the ability of ORF3b to suppress interferon induction. Thus, our findings not only help to explain the poor interferon response in COVID-19 patients but also describe the emergence of natural SARS-CoV-2 quasispecies with an extended ORF3b gene that may potentially affect COVID-19 pathogenesis.

Genome evolution of SARS-CoV-2 and its virological characteristics.

Coronavirus disease of 2019 (COVID-19), which originated in China in 2019, shows mild cold and pneumonia symptoms that can occasionally worsen and result in deaths. SARS-CoV-2 was reported to be the causative agent of the disease and was identified as being similar to SARS-CoV, a causative agent of SARS in 2003. In this review, we described the phylogeny of SARS-CoV-2, covering various related studies, in particular, focusing on viruses obtained from horseshoe bats and pangolins that belong to Sarbecovirus, a subgenus of Betacoronavirus. We also describe the virological characteristics of SARS-CoV-2 and compare them with other coronaviruses. More than 30,000 genome sequences of SARS-CoV-2 are available in the GISAID database as of May 28, 2020. Using the genome sequence data of closely related viruses, the genomic characteristics and evolution of SARS-CoV-2 were extensively studied. However, given the global prevalence of COVID-19 and the large number of associated deaths, further computational and experimental virological analyses are required to fully characterize SARS-CoV-2.